Giant spin rotation under quasiparticle-photoelectron conversion: 
Joint effect of sublattice interference and spin-orbit coupling 
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Spin- and angular-resolved photoemission spectroscopy is a basic experimental tool for unveiling 
spin polarization of electron eigenstates in crystals. We prove, by using spin-orbit coupled graphene 
as a model, that photo conversion of a quasiparticle inside a crystal into a photoelectron can be ac- 
companied with a dramatic change in its spin polarization, up to a total spin flip. This phenomenon 
is typical of quasiparticles residing away from the Brillouin zone center and described by higher rank 
spinors, and results in exotic patterns in the angular distribution of photoelectrons. 



PACS numbers: 71.70.Ej, 73.61.-r, 73.61.Wp, 79.60.-i 

Establishing spin polarization of quasiparticles in crys- 
tals is of crucial importance for semiconductor spin- 
tronics, physics of metallic surfaces, and the emerging 
field of topological insulators. The latter typically com- 
prises narrow-gap systems with strong spin-orbit cou- 
pling. Graphene, a zero-gap conductor, attracts atten- 
tion due to its linear Dirac-Weyl energy spectrum and 
prospects for applications. Its quasirelativistic spec- 
trum, manifesting itself in unconventional quantum Hall 
effect^' and Klein tunneling^, is essentially a nonrela- 
tivistic phenomenon originating from two equivalent sub- 
lattices, A and B] their effect is conveniently accounted 
for by pseudospin. However, the relativistic effects that 
entangle their spin and pseudospin degrees of freedonP 
lift the spin degeneracy of the energy spectrum. Intrinsic 
spin-orbit inte raction in graphene^ is weak, not exceed- 
ing tens of /ie V^ I ^ I ^ I^, but breaking the up-down symme- 
try by a substrate can result in a substantial extrinsic 
spin-orbit coupling. Such a coupling of the scale of 10 
meV was discovered by Varykhalov et al^ by spin- and 
angular-resolved photoemission spectroscopy (SARPES) 
techniques. Despite its modest magnitude, this coupling 
modifies essentially the zero-gap nonrelativistic spec- 
trum. This makes spin-orbit coupled graphene an excel- 
lent platform for unveiling nontrivial effects of spin-orbit 
coupling on SARPES spectra. 

In graphene the zeros of the gap are achieved in two 
nonequivalent corners of the Brillouin zone (BZ), K and 
K' [Fig. Ic]. In the vicinity of these points the non- 
relativistic quasiparticles (electrons and holes) possess a 
quasirelativistic energy spectrum E{k) = ±7/^, k be- 
ing a q uasimomentum counted relative to the K{K') 
pointPJ. Near the K point, the nonrelativistic Hamilto- 
nian is = 7((T • k), where the pseudospin cr = (<Ji, CF2) 
is a vector of Pauli matrices acting on (A, B) sublattices 
[Fig. 1(a)]. The leading term in the extrinsic spin-orbit 
coupling is TYso = ^A(cr x s);^, where s are Pauli ma- 
trices of the real spin and A is the spin-orbit coupling 
constant"^ it couples spin to the pseudospin. We dis- 
regard the breaking of the (A, B) symmetry by the stag- 
gering potential of the substrate; for some substrates it 




Figure 1: Energy spectrum and spin polarization of quasi- 
particles in spin-orbit coupled graphene. (a) Schematic of a 
graphene flake with substrate and a SARPES detector (sensi- 
tive to energy, momentum, and spin of photoelectrons). Sym- 
metry between A and B sublattices is preserved, (b) Energy 
spectrum and spin polarization near the K point. Spins are 
in-plane polarized transverse to the momentum k. The mag- 
nitude of the spin polarization of quasiparticles (arrows) is 
proportional to the group velocity. Its sign is shown for A > 0; 
it is opposite for A < 0. (c) Brillouin zone and coordinate sys- 
tems in the momentum space. 



is weak and graphene behaves as quasifreestandin^^^l^. 
The two-sublattice structure of graphene is known to re- 
sult in an interference effect in the nonrelativistic ARPES 
spectrum that produces strong photoemission anisotropy 
but does not distort the shape of the isoenergy surfaces^. 
In relativistic spectra, the joint effect of interference and 
the spin-pseudospin entanglement produces giant rota- 
tions of electron spins during photoemission resulting in 
drastic differences in the spin polarization of quasipar- 
ticles inside the crystal and photoelectrons in vacuum. 
This phenomenon is the focus of this paper. 

The Hamiltonian = TLq + TLso acts in the space 
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of four-spinors defined in tiie basis {\A ^),\A i)^\B t 
),\B I)}, tiie products of tiie two-fold degenerate Blocii 
functions of tiie K point at A{B) sublattices and spin- 
up (down) spinors. It is convenient to ciiange to a basis 
{\A T),|5 i),\B T), |A I)} and express eacii bispinor, 
^^y^, in terms of two spinors, (puj^ and 0^,^, as ^^^^ = 

I 1 . details see Ref. |15l The eigenvalues are 

E,,{k) = f (^y'J^Ti^ - n\) , (1) 

where z^, /i = ±1. The spectrum consists of two ungapped 
and two gapped hyperbolas shifted by A that are shown in 
Fig. 1(b) with their quantum numbers. The spectrum is 
similar to unbiased bilayer graphene^^, but with different 
nature of eigenstates and narrower gap^ In the new basis, 
matrices of the quasiparticle spin are S = ais, and their 
mean values in the {u^ ji) eigenstates are 

z being a unit vector perpendicular to the substrate. We 
note in passing that (S) is proportional to the group ve- 
locity dEjy^/hd'k. The sign of the chirality of eigenspinors 
{(Pun^<l^un) is determined by the sign of spin polariza- 
tion by /i, and the sign of the energy by the product u/j.. 
The polarization of all branches is depicted in Fig. 1(b). 
Such a spin-polarized spectrum with |A| ^ 13 meV was 
deduced from SARPES data taken from graphene on a 
Au/Ni(lll) substrate^^ for \Ei,^{k)\ > |A|; data for lesser 
energies are not available. 

According to Eq. ([2|, quasiparticles are in-plane po- 
larized perpendicular to the momentum k, and the mag- 
nitudes of their spins depend on k and can even van- 
ish. This is well compatible with the fact that = 3, 
because spin-orbit coupling results in large nondiagonal 
components of S in the 2x2 spin-pseudospin space, hence 
restrictions on its diagonal components ("crystal" spin) 
are relaxed. The scale of the anticipated effect follows 
from the ratio of 7io and Hso that is about k/kx^ where 
kx = |A|/27 is a spin-orbit momentum. The term TYo, 
Zeeman energy of the pseudospin in an effective mag- 
netic field k, describes the k-dependence of the interfer- 
ence of Bloch waves scattered by two equivalent sublat- 
tices, and this mechanism affects the spin sector through 
the term Hso- This spin renormalization is strongest for 
k/kx ^ 1 [Fig. 1(b)] and leads to vanishing spin at the 
symmetry points K and K' as a result of lattice inter- 
ference. We emphasize that this is impossible for simple 
lattices, where (s)^ = 1 is maintained either by in-plane 
polarizatioiP^or by out of plane spin rotation as recently 
reporte 

While the effect of sublattice interference on the an- 
gular dependence of ARPES, unveiled by S hirley et 
a/.l^, is w ell documented for monolaye^^^ and bi- 
layer CnESI graphene and for graphitJ^^'^^'^S'^, the re- 
cent progress in studying quasiparticles by SARPES 
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Figure 2: Spin-polarization of quasiparticles and photoelec- 
trons near the K point, (a) Direction (arrows) and magnitude 
(lengths) of quasiparticle spins of the lower spectrum branch; 
spins vanish for k ^ 0. (b) Photoemission probability from 
non-relativistic graphene (A = 0). Color indicates the pho- 
tocurrent /+i,_i(^) (dark = 0, bright = 2). (c) Spin polar- 
ized photocurrent from the lower spectrum branch; arrows 
indicate photoelectron spins. Polarization persists at /c ^ 
and changes fast near the "dark corridor" {k\\ > 0, k± = 0). 
For /c < 2, in units of kx. 



techniques, the novelty of results, and the large mag- 
nitude of spin-orbit coupling fo und in a number of 
system d^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^ l ^^l^^ l ^a Sol^ all call for proper un- 
derstanding of interference-generated spin patterns in 
photoemission. They are specific for the single-step 
transformation of Bloch spinors of quasiparticles into 
photoelectron plane waves with large in-plane (Brillouin) 
momentum and manifest themselves in singular k depen- 
dencies in the k ^ limit. The Fano effect in atomic 
photoionization^^ and polarization-dependent interband 
transition never show such singularities. They also are 
absent in the single-step normal photoemission showing 
remarkable polarization properties depending on the in- 
cidence conditionP^. 

For calculating photoemission, one needs to employ 
detailed wave functions that are products of the com- 
ponents of envelope spinors ^j,^ and Bloch basis func- 
tions. At the K point of the BZ [Fig. 1(c)], the \B) 
functions differ from the symmetrically equivalent to 
them \A) functions by the phase factor (— c^;), with uj = 
exp (z27r/3)^^ The uj factor reflects the effect of sub- 
lattice interference on photoemission in terms of the en- 
velope functions, while the specific form of 1^4) is of no 
importance as long as small spin-orbit corrections de- 
pending on atomic form-factors^ factorize out and can 
be disregarded. For nonrelativistic electrons (A = 0) the 
interference factor in the photoemission intensity reduces 
to [{k-/k){E/\E\) -uj] near the K point, k <^ K; here 
k- = kx — iky. The first term in the brackets originates 
from TYo, and the second from the interference of out- 
going photoelectrons (exp [iK • (Ra — Rb)] = with 
IIa,b for sublattice coordinates). The interference fac- 
tor is related to the quasiparticle pseudospin whose mean 
value equals {cix — i(Jy) = —u;{E/\E\) for a quasiparticle 
with a momentum k along K. The resulting photoemis- 
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sion is highly anisotropic as described by the large k limit 
of Eq. ([5| and displayed in Fig. 2(b). 




(arb. units) 

n-2 



Figure 3: Intensity (colorscale) and spin-polarization (ar- 
rows) for photoelectrons emitted from all four branches near 
the K point, for k < 5 and A > 0; (a) and (d) for gapped, and 
(b) and (c) for ungapped branches. Spins of photoelectrons 
and quasiparticles become identical only for /c ^ 1 and away 
from the "dark corridor" . k in units of kx . 

For relativistic electrons (A 7^ 0) the envelope spinor 
$jy^(k), describing the flux of spin-polarized photoelec- 
trons from the (z^, ja) eigenstate, can be found by adding 
the components of ^^^^(k) corresponding to the same 
spin polarizations of the basis spinors, supplied with the 
proper phase factors originating from \A) and \B) func- 
tions 



^i./i(k) oc 



(^i;](k)-c^4M(k) 

4'/](k)-c^(^l'j(k) 



(3) 



where superscripts indicate the upper and lower compo- 
nents of the spinors cp and ^. This spinor, defined with 
the accuracy to a factor depending on the atomic ma- 
trix element, intensity of the source, etc., describes the 
dependence of the spin flux on the azimuth [Fig. 1(c)]. 

The final form of ^jy^Ck) is determined by the explicit 
form of the quasiparticle spinors 



^^/i(k) = 
</>iv/i(k) = 



7/c 



1 

1 r 



(4) 



normalized as (^iy^(k)|^j,^(k)) = 1. In particular, the 
total flux equals h^{k,e) = (<l>^^(k) |<l>^^(k)) 



h^{k,e) = l + 2z//i7fc cos 19/ VA2+ 472^2^ (5) 

In the large k limit, k ^ kx^ \i coincides with the well 
known result^^ 1 + ujacosO^ with its strong depen- 
dence originating from the sublattice interference. On 



the contrary, in the small k region, /c <C /c^, emission is 
isotropic due to a strong spin-pseudospin entanglement. 

Our focus is on photoelectron spins. They are in-plane 
polarized, and in a reference system {k\\^k±) related to 
the K point [Fig. 1(c)]: 



i^{k,0) = 2fijk 



1 ^ -fk cos 0/E^^{k) 



^JW^^^fW + 2v^^k cos d 



smt 



(6) 



(1 + 7^ cosO/E^^) cos e - v\/2^k 
V'A^ +472^:^ + 2vii^k cos 



(7) 

Figures 2(a) and 2(c) show unidirectional flow for 
photoelectrons (left to right) while for quasiparticles 

Alternatively, spins can be specified by their transverse 
and longitudinal components, and s^, defined as mean 
values of (s x k)^ and (s-k), respectively; k = k//c. Then 



.(k) = 



vfiX sin^ 



^/>?^^A^^fW + 2uii'^k cos 9 



(8) 



and (sf,^)^ + ('^t/i)^ = 1- Spins of photoelectrons acquire 
a considerable longitudinal component (unless ^ = 0, tt) 
despite the fact that quasiparticles are transverse polar- 
ized. It decreases when k/k\ 00, but retains a consid- 
erable value in the moderate k/kx region. The magnitude 
of remains large even in the region of high brightness 
(6 ^ 7r/2). While for quasiparticles all spins vanish in 
the k ^ limit, for photoelectrons they do not, with 
sl^{k = 0) = z^/iAsin6>/|A| and s^^(/c = 0) = ucosO. The 
angular distribution of the total flux, depends on 

A [Figs. 2(b) and 2(c)]. 

Spins of photoelectrons coming from different branches 
are compared in Fig. 3(a-d) in a wider k range. For 
k < kx^ the difference between gapped and ungapped 
branches strikes the eye: for the former spins follow a 
simple one-sided flow, while for the latter their direction 
strongly depends on k and 6. For k > kx^ polarization 
is transverse only inside the bright sector. While the 
transverse polarization follows from the k ^ 00 expan- 
sions of Eqs. ([6| and ([7|), the convergence is nonuniform, 
and asymptotic expansions diverge along the "dark cor- 
ridor". Near it, the signs of all are opposite for quasi- 
particles and photoelectrons [Fig. 3(a) and 3(c)]. Near 
the corridor, spins rotate at the scale of the azimuths of 
(1 + uficosO) ^ {X/2jk)'^ that grows narrower with in- 
creasing k/kx' The intensity is low and vanishes when 
/c ^ 00, but the regions of the anomalous spin polariza- 
tion are wide enough to be accessible for experiment. 

Figure 4 displays energy resolved photoemission spec- 
tra that can be anticipated. Left panels of (a) and 
(b) show Lorentz- widened (width 0.2 |A|) spectra from 
all spectrum branches in two central sections, vs k±/kx 
(with k\\ = 0) and k\\/kx (with k_\_ =0). Due to A 7^ 0, 
finite intensity persists in panel (b) even inside the "dark 
corridor" . Right panels exemplify the energy dependence 
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Figure 4: FIG. 4 Energy dependence of the total flux I±(E) 
and its spin polarization P±{E) along k±. Left panels: energy 
resolved photocurrent as a function of k± (a) and k\\ (b), 
respectively. Dashed lines indicate the traces of the right 
panels. Right panels: flux I± (black) and polarization P± 
(green). See text for details. E in units of A, k in units of kx. 

of the total photocurrent I±{E) = Ip -\- Ia (black) and 
spin polarization P±{E) = {Ip — Ia)/I± (green), where 
Ip and Ia are photocurrents polarized parallel and an- 
tiparallel to a perfect Mott detector aligned along the 
positive k_\_ direction. While the details depend on the 
choice of Lorentzians, the general patterns [such as the 
sign reversal of P{E) in Fig. 4(a)] are quite general. 

The width of the spectra originating from the lifetime 
of the holes left after the photoemission is the main ob- 
stacle for resolving spin polarized spectra. Because it 
decreases with decreasing hole energy due to the elimina- 
tion of scattering channels^^, performing measurements 
close to the emission threshold should increase the spec- 
tral resolution. From this standpoint, small A provides 
advantages. 

For in-plane symmetry preserving substrates, the re- 
strictions imposed by competing mechanisms of SO cou- 
pling are essential only within a narrow vicinity of the de- 
generacy point controlled by intrinsic SO coupling [about 
10 /ieV- 10-^A (Refs. I6I7I8I9,) ] and trigonal SO correc- 



tions [for k < akj^ <^ a being a lattice constant, 
akx ~ 10~^ (Ref. llll)]; both for the parameter values of 
Ref-Hni 

It is seen from comparing Figs. 1(b) and 3(a-d) that the 
difference between the spins of quasiparticles and photo- 
electrons ((S)iy^(k) and Sjy^{k^6)) differs between three 
27r/3 sectors around the K point; for panels (a) and (b) 
of Fig. 3 it is minimal inside the internal sector that is 
bright. These data allowed reconstructing the general 
pattern of spin polarization in the k ^ kx range^^. Un- 
veiling the spin polarization in the k ^ kx range is more 
challenging. Actually, it is an open question whether re- 
constructing (S)^^(k), related to four-component spinors 
^iy^(k), from Siyjj^{k, 6>), related to two-component spinors 
<l>jy^(k), is a well posed mathematical problem, and how 
much spin information is lost during photoemission. In 
any case, reconstruction of the k < kx region should es- 
sentially include the data from external sectors. 

All calculations were performed for spin-orbit coupled 
graphene. However, the qualitative conclusions are gen- 
eral and are based on (i) a narrow-gap spectrum of quasi- 
particles with strong spin-orbit coupling described by 
four-component spinors, and the presence of the co phase 
factors that requires (ii) existence of (nearly) equivalent 
sublattices and (iii) residence of quasiparticles away from 
the BZ center. These requirements are fulfilled for the 
Bii-a.Sb^'^^ and Pbi_^Sn^Te^^ compounds. There- 
fore, we expect that their SARPES spectra of bulk elec- 
trons and the electrons from their Tamm-Shockley sur- 
face bands should manifest similar anomalies. 

In conclusion, we demonstrated a dramatic effect of the 
spin-pseudospin entanglement on the SARPES spectra of 
spin-orbit coupled graphene. Similar effects are predicted 
for different narrow-gap spin-orbit coupled materials. 
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